Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate.

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.

Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate.

In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.

Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates.

In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations.

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log10 (J/(m3 s)) = -7(5).

Intestinal variation of serotonin, melatonin, and digestive enzymes activities along food passage time through GIT in Salmo salar fed with supplemented diets with tryptophan and melatonin.

In teleosts, peripheral serotonin (5-HT) and melatonin (MEL) are synthesised in the gastrointestinal tract (GIT) and regulate secretion and motility processes. Their production is regulated by diet and the passage of food through the GIT. This study aimed to evaluate how intestinal 5-HT, melatonin, and the activity of digestive enzymes varied with food passage time through GIT in Atlantic salmon (Salmo salar). We fed fish diets supplemented with tryptophan and melatonin (L-Trp 2.5% and MEL 0.01%) and measured the activity of digestive enzymes (amylase, lipase, and total protease) in the pyloric caeca, midgut, and hindgut at different times after feeding. 5-HT levels increased in all GIT portions and diets at 120 min post-intake and were highest in the pyloric caeca. Intestinal enzymatic activity was varied with diet, post-intake time and in different intestinal portions. In conclusion, food passage time directly affects GIT 5-HT secretion and digestive enzyme activity in S. salar, and diet composition regulates S. salar GIT function.

Glucosensing capacity of rainbow trout telencephalon.

To assess the hypothesis of glucosensing systems present in fish telencephalon, we first demonstrated in rainbow trout, by in situ hybridisation, the presence of glucokinase (GK). Then, we assessed the response of glucosensing markers in rainbow trout telencephalon 6 hours after i.c.v. treatment with glucose or 2-deoxyglucose (inducing glucoprivation). We evaluated the response of parameters related to the mechanisms dependent on GK, liver X receptor (LXR), mitochondrial activity, sweet taste receptor and sodium-glucose linked transporter 1 (SGLT-1). We also assessed mRNA abundance of neuropeptides involved in the metabolic control of food intake (agouti-related protein, neuropeptide Y, pro-opiomelanocortin, and cocaine- and amphetamine-related transcript), as well as the abundance and phosphorylation status of proteins possibly involved in linking glucosensing with neuropeptide expression, such as protein kinase B (AkT), AMP-activated protein kinase (AMPK), mechanistic target of rapamycin and cAMP response element-binding protein (CREB). The responses obtained support the presence in the telencephalon of a glucosensing mechanism based on GK and maybe one based on LXR, although they do not support the presence of mechanisms dependent on mitochondrial activity and SGLT-1. The mechanism based on sweet taste receptor responded to glucose but in a converse way to that characterised previously in the hypothalamus. In general, systems responded only to glucose but not to glucoprivation. Neuropeptides did not respond to glucose or glucoprivation. By contrast, the presence of glucose activates Akt and inhibits AMPK, CREB and forkhead box01. This is the first study in any vertebrate species in which the response to glucose of putative glucosensing mechanisms is demonstrated in the telencephalon. Their role might relate to processes other than homeostatic control of food intake, such as the hedonic and reward system.

Computational study of the interplay between intermolecular interactions and CO2 orientations in type I hydrates.

Carbon dioxide (CO2) molecules show a rich orientation landscape when they are enclathrated in type I hydrates. Previous studies have described experimentally their preferential orientations, and some theoretical works have explained, but only partially, these experimental results. In the present paper, we use classical molecular dynamics and electronic density functional theory to advance in the theoretical description of CO2 orientations within type I hydrates. Our results are fully compatible with those previously reported, both theoretical and experimental, the geometric shape of the cavities in hydrate being, and therefore, the steric constraints, responsible for some (but not all) preferential angles. In addition, our calculations also show that guest-guest interactions in neighbouring cages are a key factor to explain the remaining experimental angles. Besides the implication concerning equation of state hydrate modeling approximations, the conclusion is that these guest-guest interactions should not be neglected, contrary to the usual practice.

On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation.

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

Understanding the Phase Behavior of Tetrahydrofuran + Carbon Dioxide, + Methane, and + Water Binary Mixtures from the SAFT-VR Approach.

The high-pressure phase diagrams of the tetrahydrofuran(1) + carbon dioxide(2), + methane(2), and + water(2) mixtures are examined using the SAFT-VR approach. Carbon dioxide molecule is modeled as two spherical segments tangentially bonded, water is modeled as a spherical segment with four associating sites to represent the hydrogen bonding, methane is represented as an isolated sphere, and tetrahydrofuran is represented as a chain of m tangentially bonded spherical segments. Dispersive interactions are modeled using the square-well intermolecular potential. In addition, two different molecular model mixtures are developed to take into account the subtle balance between water-tetrahydrofuran hydrogen-bonding interactions. The polar and quadrupolar interactions present in water, tetrahydrofuran, and carbon dioxide are treated in an effective way via square-well potentials of variable range. The optimized intermolecular parameters are taken from the works of Giner et al. (Fluid Phase Equil. 2007, 255, 200), Galindo and Blas (J. Phys. Chem. B 2002, 106, 4503), Patel et al. (Ind. Eng. Chem. Res. 2003, 42, 3809), and Clark et al. (Mol. Phys. 2006, 104, 3561) for tetrahydrofuran, carbon dioxide, methane, and water, respectively. The phase diagrams of the binary mixtures exhibit different types of phase behavior according to the classification of van Konynenburg and Scott, ranging from types I, III, and VI phase behavior for the tetrahydrofuran(1) + carbon dioxide(2), + methane(2), and + water(2) binary mixtures, respectively. This last type is characterized by the presence of a Bancroft point, positive azeotropy, and the so-called closed-loop curves that represent regions of liquid-liquid immiscibility in the phase diagram. The system exhibits lower critical solution temperatures (LCSTs), which denote the lower limit of immiscibility together with upper critical solution temperatures (UCSTs). This behavior is explained in terms of competition between the incompatibility with the alkyl parts of the tetrahydrofuran ring chain and the hydrogen bonding between water and the ether group. A minimum number of unlike interaction parameters are fitted to give the optimal representation of the most representative features of the binary phase diagrams. In the particular case of tetrahydrofuran(1) + water(2), two sets of intermolecular potential model parameters are proposed to describe accurately either the hypercritical point associated with the closed-loop liquid-liquid immiscibility region or the location of the mixture lower- and upper-critical end-points. The theory is not only able to predict the type of phase behavior of each mixture, but also provides a reasonably good description of the global phase behavior whenever experimental data are available.

Effects of insulin treatment on the response to oleate and octanoate of food intake and fatty acid-sensing systems in rainbow trout.

We hypothesized that food intake and the response of fatty acid (FA)-sensing systems in hypothalamus, liver, and Brockmann bodies of rainbow trout to raised levels of oleate (OL) or octanoate (OCT) is modified by insulin treatment. To assess this hypothesis, 15 fish per group received intraperitoneally 10-mL/kg injection of saline solution alone (control), or containing insulin (2-mg bovine insulin/kg body mass), OL (300 µg/kg), OCT (300 µg/kg), insulin + OL, or insulin + OCT to be sampled 6 h later to assess parameters related to FA sensing. Our results suggest that the modulatory role of insulin on the responses of hypothalamic FA-sensing systems to changes in circulating levels of OL or OCT was of minor importance in contrast to the mammalian model. However, this is in contrast with the effects observed in another experiment assessing changes in food intake after similar treatments because insulin treatment enhanced the anorectic effects of FA alone, and the effect was especially relevant (P < 0.001) for OCT, in contrast with the mammalian model where this FA is not inducing an anorectic response. In liver and Brockmann bodies, insulin treatment enhanced the responses to OL or OCT treatment in parameters related to FA sensing. Therefore, we provide for the first time in fish, and in a non-mammalian vertebrate, evidence for the modulation of FA-sensing systems by insulin.

Molecular dynamics simulation of CO2 hydrates: Prediction of three phase coexistence line.

The three phase equilibrium line (hydrate-liquid water-liquid carbon dioxide) has been estimated for the water + carbon dioxide binary mixture using molecular dynamics simulation and the direct coexistence technique. Both molecules have been represented using rigid nonpolarizable models. TIP4P/2005 and TIP4P/Ice were used for the case of water, while carbon dioxide was considered as a three center linear molecule with the parameterizations of MSM, EPM2, TraPPE, and ZD. The influence of the initial guest occupancy fraction on the hydrate stability has been analyzed first in order to determine the optimal starting configuration for the simulations, paying attention to the influence of the two different cells existing in the sI hydrate structure. The three phase coexistence temperature was then determined for a pressure range from 2 to 500 MPa. The qualitative shape of the equilibrium curve estimated is correct, including the high pressure temperature maximum that determines the hydrate re-entrant behaviour. However, in order to obtain quantitative agreement with experimental results, a positive deviation from the classical Lorentz-Berthelot combining rules must be considered.

Gradation of the stress response in rainbow trout exposed to stressors of different severity: the role of brain serotonergic and dopaminergic systems.

After an intense acute stressor, fish develop a metabolic and behavioural response that usually lasts for several hours. Brain monoaminergic systems, particularly the serotonergic system, appear to play a key role in the central regulation of the stress response. However, the influence of stressor severity on brain monoaminergic systems and on the induced stress responses is yet poorly understood. We hypothesise that serotonergic system could have a direct role in the integration of sensory information during stressor exposure and in the organisation of the subsequent integrated stress response. According to our hypothesis, a low stressor intensity would induce a low response of brain serotonergic system and therefore stress responses of low magnitude and duration. To test this hypothesis, we exposed fish to handling disturbance for 5 s, 15 s or 3 min. We sampled fish at 0 (controls), 3, 15, 45 and 240 min after the start of the stress protocol. Brain levels of serotonin, dopamine and their respective main oxidative metabolites were quantified, along with plasma levels of stress markers (catecholamines, cortisol, glucose and lactate). Regarding stress markers, the 5-s and 15-s stress protocols induced similar and relatively low elevations in all parameters assessed. As expected, the 3-min protocol induced responses of a higher intensity and duration in all plasma parameters. Interestingly, the alterations of brain monoaminergic systems did not follow the same trend. The three stress protocols induced increases in the serotonergic activity in all brain regions analysed (hypothalamus, telencephalon and medulla oblongata), independently of the duration of the handling disturbance, whereas the effects on the dopaminergic system were minor and brain region-dependent. These data suggest that the brain serotonergic system, although likely involved in the recognition of the stressor stimuli, is not the only actor determining the magnitude and duration of the acute stress response in trout.

Oral administration of melatonin counteracts several of the effects of chronic stress in rainbow trout.

To assess a possible antistress role of melatonin in fish, we orally administered melatonin to rainbow trout for 10 d and then kept the fish under normal or high stocking density conditions during the last 4 d. Food intake; biochemical parameters in plasma (cortisol, glucose, and lactate concentrations); liver (glucose and glycogen concentrations, and glycogen synthase activity); enzyme activities of amylase, lipase, and protease in foregut and midgut; and content of the hypothalamic neurotransmitters dopamine and serotonin, as well as their oxidized metabolites, 3,4-dihydroxyphenylacetic acid and 5-hydroxy-3-indoleacetic acid, were evaluated under those conditions. High stocking density conditions alone induced changes indicative of stress conditions in plasma cortisol concentrations, liver glycogenolytic potential, the activities of some digestive enzymes, and the 3,4-dihydroxyphenylacetic acid-to-dopamine and 5-hydroxy-3-indoleacetic acid-to-serotonin ratios in the hypothalamus. Melatonin treatment in nonstressed fish induced an increase in liver glycogenolytic potential, increased the activity of some digestive enzymes, and enhanced serotoninergic and dopaminergic metabolism in hypothalamus. The presence of melatonin in stressed fish resulted in a significant interaction with cortisol concentrations in plasma, glycogen content, and glycogen synthase activity in liver and dopaminergic and serotoninergic metabolism in the hypothalamus. In general, the presence of melatonin mitigated several of the effects induced by stress, supporting an antistress role for melatonin in rainbow trout.

Arginine vasotocin treatment induces a stress response and exerts a potent anorexigenic effect in rainbow trout, Oncorhynchus mykiss.

The peptide arginine vasotocin (AVT), homologous to mammalian arginine vasopressin, is involved in many aspects of fish physiology, such as osmoregulation, regulation of biological rhythms, reproduction, metabolism or responses to stress, and the modulation of social behaviours. Because a decrease in appetite is a general response to stress in fish and other vertebrates, we investigated the role of AVT as a possible food intake regulator in fish. We used i.c.v. injections for central administration of AVT to rainbow trout (Oncorhynchus mykiss). In a first experiment, we evaluated the temporal response of food intake after AVT treatment. In a second experiment, we investigated the effects of central AVT administration on the response of typical stress markers (plasma cortisol, glucose and lactate), as well as brain serotonergic, noradrenergic and dopaminergic activity. In addition, the mRNA levels of genes involved in food intake regulation [neuropetide Y, pro-opiomelanocortin (POMC), cocaine- and amphetamine-regulated transcript (CART) and corticotrophin-releasing factor (CRF)] and in CRF- (CRF-binding protein) and AVT-signalling (pro-VT and AVT receptor), were also assessed after AVT treatment. Our results showed that AVT is a potent anorexigenic factor in fish. Increases of plasma cortisol and glucose after AVT treatment strongly suggest that AVT administration induced a stress response and that AVT action was mediated by hypothalamic-pituitary-interrenal axis activation, which was also supported by the increase of the serotonergic activity in trout telencephalon and hypothalamus. The increased hypothalamic levels of POMC and CART suggest that these peptides might have a role in the anorexigenic action of AVT, whereas the involvement of CRF signalling is unclear.

Influence of the long-range corrections on the interfacial properties of molecular models using Monte Carlo simulation.

We analyze the influence of the long-range corrections, due to the dispersive term of the intermolecular potential energy, on the surface tension using direct simulation of the vapour-liquid interface of different molecular models. Although several calculation methods have been proposed recently to compute the fluid-fluid interfacial properties, the truncation of the intermolecular potential or the use of the tail corrections represents a contribution relevant from a quantitative perspective. In this work, a simplified model for methane, namely a spherical Lennard-Jones intermolecular potential, has been considered first, and afterwards other models including rigid non polarizable structures with both Lennard-Jones sites and point electric charges, representing some of the most popular models to describe water (namely the original TIP4P model, and the TIP4P/Ew and TIP4P/2005 versions), and carbon dioxide (MSM, EPM2, TraPPE, and ZD models) have been studied. Our results show that for all cases tested, including those in which the electrostatic interactions may be predominant, an incomplete account of the long-range corrections produces a systematic underestimation of the computed interfacial tension.

Vapor-liquid interfacial properties of rigid-linear Lennard-Jones chains.

We have obtained the interfacial properties of short rigid-linear chains formed from tangentially bonded Lennard-Jones monomeric units from direct simulation of the vapour-liquid interface. The full long-range tails of the potential are accounted for by means of an improved version of the inhomogeneous long-range corrections of Janecek [J. Phys. Chem. B 110, 6264-6269 (2006)] proposed recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] valid for spherical as well as for rigid and flexible molecular systems. Three different model systems comprising of 3, 4, and 5 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test-area method. In addition to the surface tension, we also obtain density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length and rigidity on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The surface tension has been scaled by critical properties and represented as a function of the difference between coexistence densities relative to the critical density.

On interfacial tension calculation from the test-area methodology in the grand canonical ensemble.

We propose the extension of the test-area methodology, originally proposed to evaluate the surface tension of planar fluid-fluid interfaces along a computer simulation in the canonical ensemble, to deal with the solid-fluid interfacial tension of systems adsorbed on slitlike pores using the grand canonical ensemble. In order to check the adequacy of the proposed extension, we apply the method for determining the density profiles and interfacial tension of spherical molecules adsorbed in slitlike pore with different pore sizes and solid-fluid dispersive energy parameters along the same simulation. We also calculate the solid-fluid interfacial tension using the original test-area method in the canonical ensemble. Agreement between the results obtained from both methods indicate that both methods are fully equivalent. The advantage of the new methodology is that allows to calculate simultaneously the density profiles and the amount of molecules adsorbed onto a slitlike pore, as well as the solid-fluid interfacial tension. This ensures that the chemical potential at which all properties are evaluated during the simulation is exactly the same since simulations can be performed in the grand canonical ensemble, mimicking the conditions at which the adsorption experiments are most usually carried out in the laboratory.

Ghrelin effects on central glucosensing and energy homeostasis-related peptides in rainbow trout.

Although the role of ghrelin (GHRL) on fish appetite regulation had been widely studied in past years, its involvement in the regulation of glucose metabolism had been little explored. In the present study we hypothesize that GHRL may have a role in the regulation of glucose homeostasis in fish. Therefore, we carried out different experimental approaches in rainbow trout to assess brain glucosensing potential and glucose metabolism in response to GHRL treatment. We found that after either systemic or central GHRL administration to trout deprived of food, glycemia remained unaffected, whereas (in clear contrast with the mammalian model) a consistent activation of the main glucosensing markers (glucose transporter 2, glucokinase, and ATP-sensitive inward rectified K+ channel) was noticed in both hypothalamus and hindbrain. Some of these results were further confirmed by in vitro incubations of hypothalamus and hindbrain in the presence of GHRL. Despite the lack of changes in glycemia, we suggest that the changes elicited by GHRL on the glucosensing system are direct and could be related to a helper action of this hormone when glucose arrived in the postprandial phase. Moreover, we also studied the effect of GHRL treatment on the expression of several food intake-related neuropeptides, such as neuropeptide Y, cocaine- and amphetamine-regulated transcript, pro-opiomelanocortin, and corticotropin-releasing factor. We observed an important variability in the effects of GHRL attributable either to the route of GHRL administration or to the brain regions assessed, which could help explain the contradictory results described in fish literature about GHRL role in food intake control.

An examination of the ternary methane + carbon dioxide + water phase diagram using the SAFT-VR approach.

In this work, the molecular based Variable Range Statistical Associating Fluid Theory (SAFT-VR) has been used to estimate the global phase equilibria diagram of the ternary mixture water + carbon dioxide + methane, over a wide pressure and temperature range. An accurate determination of the phase equilibria of this mixture is relevant in Petrophysics, as, for instance, in enhanced natural gas recovery from low permeability reservoirs (the so-called tight gas reservoirs), or in geology, as it is the basic composition of many geological fluids. A previous study on the phase behavior of the binary mixtures involved is presented, using in a transferable manner the characteristic molecular parameters for the three molecules involved. The ternary mixture presents a very rich and complex phase behavior, with a wide region of the thermodynamic space of phases (at higher pressures) presenting a large gap of ternary liquid-liquid equilibria, that upon descending pressures leads to the transition to a three-phase liquid-liquid-vapor equilibria region, and both regions are separated by a continuous critical end point line. The ability of the theory to describe this complex multicomponent mixture phase transition with a reduced and physically sound set of characteristic parameters must be underlined.

Evidence for a gut-brain axis used by glucagon-like peptide-1 to elicit hyperglycaemia in fish.

In mammals, glucagon-like peptide-1 (GLP-1) produces changes in glucose and energy homeostasis through a gut-pancreas-brain axis. In fish, the effects of GLP-1 are opposed to those described in other vertebrates, such as stimulation of hyperglycaemia and the lack of an effect of incretin. In the present study conducted in a teleost fish such as the rainbow trout, we present evidence of a gut-brain axis used by GLP-1 to exert its actions on glucose and energy homeostasis. We have assessed the effects of GLP-1 on glucose metabolism in the liver as well as the glucose-sensing potential in the hypothalamus and hindbrain. We confirm that peripheral GLP-1 administration elicits sustained hyperglycaemia, whereas, for the first time in a vertebrate species, we report that central GLP-1 treatment increases plasma glucose levels. We have observed (using capsaicin) that at least part of the action of GLP-1 on glucose homeostasis was mediated by vagal and splanchnic afferents. GLP-1 has a direct effect in parameters involved in glucose sensing in the hindbrain, whereas, in the hypothalamus, changes occurred indirectly through hyperglycaemia. Moreover, in the hindbrain, GLP-1 altered the expression of peptides involved in the control of food intake. We have elaborated a model for the actions of GLP-1 in fish in which this peptide uses a mammalian-like ancestral gut-brain axis to elicit the regulation of glucose homeostasis in different manner than the model described in mammals. Finally, it is worth noting that the hyperglycaemia induced by this peptide and the lack of incretin function could be related to the glucose intolerance observed in carnivorous teleost fish species such as the rainbow trout.